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Abstract 

We investigate the quench of Ising and Potts models via Monte Carlo dynamics, and find that 
the distribution of the site-site interaction energy has the same form as in the equilibrium case. 
This form directly derives from the Boltzmann statistics and allows to measure the instantaneous 
temperature during the systems relaxation. We find that, after an undercritical quench, the system 
equilibrates in a finite time at the heatbath temperature, while the energy still decreases due to 
the coarsening process. 
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I. INTRODUCTION 



A key quantity in Monte Carlo simulation of equilibrium systems is the heathbath tem- 
perature T. Monte Carlo dynamics follows a Markovian process whose stationary state is 
proportional to the Boltzmann probability density at the chosen temperature, i.e.: 

P(E) ocexp[-/?£], 

where E is the system energy and f3 is the inverse heathbath temperature in Boltzmann's 
constant units. Temperature is usually a parameter with an imposed value, the fluctuating 



quantity being t 
recently shown 



;he energy, and this corresponds to simulate a canonical ensemble. It has been 
"\ that also a microcanonical formulation can be constructed, in which the 
system energy is constrained whereas the temperature is computed as a dependent quantity. 
In this formulation the temperature is found to assume the correct value at equilibrium, thus 
establishing the equivalence of simulating the microcanonical and the canonical ensemble. 

Our work aims to investigate the possibility of associating an istantaneous temperature 
to a system even when this is out of the quilibrium. Therefore we have employed the method 
developed in Ref. Q in a non equilibrium case, by computing the temperature of the Ising 
model along a quench from high temperature. We have then extended the method to the 
Potts model with q states. With similar aims different approaches for off-lattice systems 
(molecular dynamics) have been recently introduced that however are not suitable 

for lattice models. 

In the method adopted here computation bases on the statistical distribution of the site- 
site interaction energies. The most interesting finding is that the shape of the distribution 
computed along the quench is very similar to the one which characterizes the system at 
equilibrium, but corresponds to a different value of the temperature. Thence it is possible 
to associate an instantaneous temperature to the non equilibrium states of the system. 

In Section II we recall how temperature can be determined in the Ising system at equi- 
librium; in Sec. Ill the method is applied to the same system in non equilibrium and in Sec. 
IV it is extended to the Potts model with q states. Main results are summarized in Sec. V. 
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II. THE ISING MODEL 

In Ref. |l| the equilibrium temperature is computed on the base of the statistical prop- 
erties of the site-site interaction energies for a generic lattice gas model. In the case of the 
Ising model, with Hamiltonian 

H = -J2 Wj, (!) 

(i<3) 

where the sum (. . .) is over nearest neighbors sites and — ±1, the temperature is obtained 
from the distribution of the quantity: 

7 = -Y. a o a v, 

M 

that represents the interaction energy of a generic site k = with its neighbors [v]. The 
Hamiltonian, Eq. (1), can be separated into H = 7 + H r , where 

7~Lr = — E a i a m, 

(Km) 

with the the sum extended to all pairs not including the site 0. Say P({a}) the (Boltzmann) 
probability for the system of N + 1 sites to stay in the configuration {a} = {a ,cr 1 , . . . ,a N }, 
this is given by 

where (3 is the inverse temperature (in Boltzmann's constant units) and Z the partition 
function. The probability -P(e) of a configuration with 7 = e, (e = —4,-2,0,2,4) can be 
obtained as 

P(l = e)=E 7)^(M) = E 6 ^ 7) expb3 7 ] exp[-/^ r ] = 
M M 

exp[-/9e] E J2 5 ( e '7) exp[-[3H r }, 

a {a r } 

where 5(a, (3) — 1 if a — (3 and = otherwise, and {a r } is the set of all possible configurations 
of the system with o"o fixed. All and only those configurations yielding 7 = e contribute to the 
sums in the last term, so they add to a number which only may depend on the temperature 
and e, thus: 

P(7 = e) = e-*a09,e). 



3 



In order to avoid the explicit computation of a(/3, e), which is as difficult as that of the 
partition function, let us consider: 

P{l = -e)=E 5 (- e > = E <K- e > 7) exp[-/3 7 ] exp[-/3W r ] = 

M M 

ex p[/^ e ] E E 5 ( _e ' 7) exp[-pn r ]. 

°"0 {oy} 

By letting cr — > er = — (To, thus 7 — > —7, and the last expression yields 
P(l = -e) = exp[/3e] E E <K- e ' "7) exp[-/W r ], 

O"0 {°"r} 

that is, since 5(— e, —7) = <5(e, 7), 

P( 7 = _ e )=e^a(/3,e). 

Thence the ratio 

P (7 = g) = 



P( 7 = -e) 
does not depend on a(e, 

This equation derives from the canonical distribution and can be used to estimate the 
lattice temperature Tg in numerical simulations of the Ising model since, by the law of the 
large numbers, one can approximate 

P( 7 = -e) " n( 7 = -e) 1 J 

being n( 7 = e) the number of lattice sites with 7 = e, and thence 

r = T B -e (3) 

with 

r = -2 (In iTT 1 . 

If necessary, better estimates can be obtained by computing averages of n(e, (3) over 
time, since the system is at equilibrium. In Ref. |l( it has been checked that Eq. (jSJ) 
is well satisfied when simulating the canonical ensemble, and has been used to compute 
temperature in a microcanonical Monte Carlo simulation with the Kawasaki 5} dynamics. 
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III. NON EQUILIBRIUM ISING MODEL 



Monte Carlo simulations are extensively employed even for studying non equilibrium 
properties of systems j^. Among the others, the properties of systems prepared far from 
equilibrium and then allowed to relax towards their equilibrium state. If for instance the 
system is initially at high temperature and is put then in contact with a finite temperature 
heat bath, it starts to cool down and eventually approaches the heath bath temperature. 
Despite the intrinsic artificial nature of Monte Carlo dynamics, physical relevance is generally 
attributed to results obtained via single site dynamics 0. It is therefore tempting to use 
Eq. (2) for probing the system temperature during the cooling. One wishes, in particular, 
to answer the following two questions: 

1) is the system at equilibrium on very short time scales? i.e. does Eq. (J3J), which was 
derived by assuming the Boltzmann (equilibrium) statistics, hold during the system 
cooling? 

2) if yes, how does temperature decrease in time during the cooling? 

In order to answers the above questions we have performed Monte Carlo simulations in 
which the 2d Ising model was quenched from infinite to finite temperature. We have found 
that the answer to question 1) is affirmative. 

In Fig. 1, r(e) is shown for a quench from T = oo to T = 3.0 at some different times. 
In this case the quench temperature is above T c , and a few Monte Carlo steps are sufficient 
for driving the system to thermal equilibrium. However, before this happens the system 
already displays a linear relation between e and r, from which an istantaneous temperature 
can be derived, as seen also in the inset where r/e is shown. A linear regression of r vs e at 
equilibrium yields T = 2.996, with correlation coefficient = 1. 

A subcritical quench is shown in Fig 2. Here T = 2.0. In this case the system takes a long 
time to reach thermal equilibrium, but Eq. (jSJ) is quite well satisfied. It can be seen from 
the inset that deviations from the linear dependence exist. They are however small for any 
e, and their average is zero. Quasi-equilibrium is satisfied at any instant and in this case a 
linear regression on the last curve gives Tg = 2.058 with a correlation coefficient = 0.9998. 

The existence of such an equilibrium-like distribution for the system allows for extracting 
a Boltzmann temperature Tg during the system coooling. Figure 3 shows the behavior of 
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the energy as function of this temperature for a quench at T = 2.0, i.e. below the critical 
temperature T c ~ 2.26. It is seen that a first regime exists in which the two quantities are 
proportional. Then, a second regime is entered where a remarkable property can be noticed: 
the Boltzmann temperature attained that of the heatbath but, on the contrary, the energy 
is still relaxing. In fact this kind of systems is characterized by a coarsenig process and 
takes an infinite time to go to equilibrium, the energy decaying as 1/^JJt), but it is seen 
that temperature is already at equilibrium during the coarsening. This behaviour is better 
pointed out in the inset, where it is seen that energy decreases algebraically while Tg is 
constant. 

When the system is quenched at a temperature above T c , energy relaxes exponentially 
to equilibrium and the second regime does not appear. This is shown in Fig. 4 and in the 
related inset. For comparison, the equilibrium curve is also reported in the figures, showing 
that it still lies below the non-equilibrium one. 

IV. THE POTTS MODEL 

In order to test the above conclusions on a wider class of systems, we have derived Eq. © 
for the case of the q states Potts model, with Hamiltonian: 



where again the sum is over the nearest neighbors, 5 a p is the Kronecker's function, and 7]k 
is the state of site k. Each site can assume one out of q different states, so that one can 
identify the state with an integer: r\ k = 1,2, ... ,q and the system energy Eq. can be 
written as 



h= £(i-<W 



(4) 



(*<?') 



E{rj) =j + H 



where 



7 = £(1 - Kvv) 



(5) 



V 



is the interaction energy of the site in with its neighbors, and 



(i<3) 



is the energy interaction of the rest of the system. 
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Let us now define 7 P the energy Eq. (J3j) when r]o = p and similarly 7^ the same quantity 
when rjo = u. In analogy with the case of the Ising model: 

— r = exp[-/3(7 P - 7^)J = exp[-/37 P ^J. (6) 

Thus the probability for 7 pa; to assume the value e is given by 

Yl 8(e,'y pu )P(j),7ii,...,Ti N ) = 

X] ^( e ; 7 P ^)exp[-/57 p jP(u;,?7i, ...,77Ar) = 

1i,-,1iv 

exp[-/3e] 51 d(t,l P u J )P(u,rii,---,VN)- 

rn,...,ijN 

By multiplying the left hand side by 1 = J2 V0 < K r /o ) P) an d he right hand side by 1 = 
E,, 5(rj Q ,u) one gets 

where {77} is the system configuration, and thus finally 

exp[-/3 £ ] = P) 
(5(r7o,w)5(e,7 pw )) 

Thence even for the Potts model the temperature can be derived from the site-site energy 

statistics. The same argument as above holds for any pair of different states rj^ , T]q , so 
that one can also write 

1 _ (^(^o,^o a V(e,7„M„(«)) 

exp [-M = Y m . 8 ) 

This can be helpful for improving statistics in out of equilibrium simulations where the 
thermal average cannot be obtained by averaging over time like in equilibrium cases. One 
thus replaces thermal averages with system averages in Eq. (jHJ). 

It is found that the Potts model obeys Eq. (refpottstemperature) aven in the non equi- 
librium regime, allowing to extract an instantaneous temperature Tg like done for the Ising 
model. Figure 5 shows the Boltzmann temperature Tb vs time for a quench of a seven state 
Potts model at T = 0.5, i.e. below the critical temperature T c ~ 0.77. It is seen that like 
for the Ising model a first regime exixts in which the two quantities are proportional. Then 
even in this case the second regime is entered in which Tb is that of the heatbath but the 
energy is still relaxing due to the coarsening process. The supercritical case reported in 
Fig. 6 shows that T B — T after a certain time, as for the Ising model. 



V. LOW TEMPERATURE BEHAVIOUR 



We have performed many quenches at different final temperature T for the the Potts 
model with q = 7 and q = 2. The latter is equivalent to the Ising model by shifting energies 
by adding 2 and compressing them by a factor 1/2 (T — c ~ 1.13). A common feature is 
that as T is decreased below a given value, the system seems unable to reach the equilibrium 
temperature. 

Figure 7 shows Tb vs time for the Ising model with linear lattice size L = 1000. The 
Boltzmann temperature attains the value of T if this is above TJ > 0.3 but if T is equal or 
below this value it doesn't. We expect that this is a finite size effect and that the system 
eventually relaxes to the ground state with Tg = 0. The case of the Potts model shown 
in Fig. 8. In this case T\ ~ 0.3. It can be observed that T\ « 0.42T C for Potts whereas 
T\ ph 0.15T C for Ising. Moreover the limit energy is clearly larger for Potts tha for Ising. 
Figure 9 displays the behaviour of Tb after a quench at T = 0. An increase of T\ with 
increasing size for both Ising and Potts model is seen. Although extrapolation to L = oo 
appears to give a still finite temperature, it must be point out that the evaluation of Tb via 
Eq. (j2J) becomes difficult for quenches at very low temperature, since almost all sites are in 
the same state and statistical sampling for e ^ becomes very poor. Reliable estimates of 
Tb can be obtained only for finite values of the quench temperature above a bound which 
depends on system size and kind. 



VI. SUMMARY 



In this work a method for associating an instantaneous temperature to non equilibrium 
lattice models has been presented which is based on the the probability distribution of the 
site-site interaction energies and which is derived from the equilibrium Boltzmann statistics. 
It is found that the same form of distribution is obeyd by the system while cooling after a 
quench, allowing to associate a Boltzmann-like temperature to its non equilibrium states. 
While cooling, the Boltzmann temperature decreases until the value of the heath bath is 
attained, both for supercritical and undercritical quenches. The latter case is particularly 
worth of notice, since there the system undergoes a coarsening process and equilibrium is 
never attained from the energetical point of view The method proofs valid for generic Potts 
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model with q states, provided the heath bath temperature is not to low and the statistical 
sampling is effective for the considered system size. 
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FIG. 1: The distribution of r(e), Eq. (j3J), for the Ising model after a supercritical quench at 
different times. In the inset the quantity r(e)/e) is shown (see text). The time interval between 
curves is of 5 mcs. 




FIG. 2: The distribution of r(e), Eq. Q, for the Ising model out of equilibrium. In the inset the 
quantity r/e is shown. The time interval between curves is of 5 mcs. 
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FIG. 3: Energy vs Temperature during a quench of the Ising model at T < T c (black curve). The 
red dashed curve is the equilibrium curve. Inset: energy and Tb vs time. 
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FIG. 4: Energy vs Temperature during a quench of the Ising model at T > T c (black curve). The 
red dashed curve is the equilibrium curve. Inset: energy and Tg vs time. 
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FIG. 5: Energy vs Temperature during a quench of the 7 state Potts model at T < T c (black 
curve). The red curve is the equilibrium curve. The inset shows and energy behavior in time. 
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FIG. 6: Energy vs Temperature during a quench of the 7 state Potts model at T > T c (black 
curve). Red curve is the equilibrium curve. In the inset, and energy vs time. 
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FIG. 9: Energy vs Boltzmann temperature for a quench of the Ising and Potts models at T = 
in lattices of different size. 
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